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Abstract. In this paper we study the effect of cosmic neutrino decoupling on 
the spectrum of cosmological gravitational waves (GWs). At temperatures T 2> 
1 MeV, neutrinos constitute a perfect fluid and do not hinder GW propagation, 
while for T <C 1 MeV they free-stream and have an effective viscosity that 
damps cosmological GWs by a constant amount. In the intermediate regime, 
corresponding to neutrino decoupling, the damping is frequency-dependent. GWs 
entering the horizon during neutrino decoupling have a frequency / ~ InHz, 
corresponding to a frequency region that will be probed by Pulsar Timing Arrays 
(PTAs). In particular, we show how neutrino decoupling induces a spectral feature 
in the spectrum of cosmological GWs just below 1 nHz. We briefly discuss the 
conditions for a detection of this feature and conclude that it is unlikely to be 
observed by PTAs. 



PACS numbers: 04.30.Nk, 95.85.Sz, 14.60.Lm 

1. Introduction 

The presence in the Universe today of a stochastic background of gravitational waves 
(GWs) is a quite general prediction of several early cosmology scenarios. In fact, the 
production of gravitational waves is the outcome of many processes that could have 
occurred in the early phases of the cosmological evolution, like the amplification of 
vacuum fluctuations in inflationary [1] and pre-big-bang cosmology scenarios [2 , phase 
transitions [3] , and finally the oscillation of cosmic strings loops [4j . In most of these 
cases, the predicted spectrum of gravitational waves extends over a very large range 
of frequencies. 

The detection of such primordial gravitational waves, produced in the early 
Universe, would be a major breakthrough in cosmology and high energy physics. This 
is because gravitational waves decouple from the cosmological plasma at very early 
times, when the temperature of the Universe is of the order of the Planck energy. In 
this way, relic gravitational waves provide us a "snapshot" of the Universe near the 
Planck time, in a similar way as the cosmic microwave background (CMB) radiation 
images the Universe at the time of recombination. 
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The extremely low frequency region (/ < 10~ 15 Hz) in the spectrum of primordial 
gravitational waves can be probed through the anisotropics of the CMB. In particular, 
gravitational waves leave a distinct imprint in the so-called magnetic or B-modes of 
its polarization field [S] [B]. The frequency region between ~ 10 Hz and few kHz 
is probed by operating large scale interferometric GW detectors like LIGO [7] and 
VIRGO [8], that, although designed with the aim to detect astrophysical signals, can 
possibly also detect signals of cosmological origin [9] . The LISA space interferometer 
[TO] , that will hopefully operate in the 2020s, will probe the frequency region between 
10~ 4 and 1 Hz. Finally, pulsar observations can be used to obtain information on 
the stochastic GW background, through the technique known as pulsar timing. The 
so-called "Pulsar Timing Arrays" will probe the region from ~1 to ~ 100 nHz, with 
a maximum sensitivity between ~3 and ~10 nHz [11] . 

In order to compare the theoretical predictions with the expected instrument 
sensitivities, one needs to evolve the GWs from the time of their production to the 
present. This is done by assuming that gravitational waves propagate freely across 
the Universe. In this case, the only effect is that the amplitude of the wave decreases 
due to the expansion of the Universe. However, GWs are sourced by the anisotropic 
stress part of the energy-momentum tensor of matter, so that the above assumption 
is well-motivated only when anisotropic stress can be neglected. It is already known 
that the anisotropic stress of free streaming neutrinos acts as an effective viscosity, 
absorbing gravitational waves in the low frequency region, thus resulting in a damping 
of the B-modes of CMB [HI |13j EH OH Q3] 

In the present work we aim to study a possible signature of neutrino decoupling 
on the spectrum of cosmological GWs. During decoupling, the effective viscosity of 
neutrinos is increasing from zero (its value at high temperatures, when the neutrinos 
are tightly coupled to the cosmological plasma) to a finite value. Neutrino decoupling 
happens when the temperature of the Universe is T ~ 1 MeV, or redshift z ~ 10 10 ; 
GWs entering the horizon at that time have a frequency / ~ 1 nHz. Since the 
viscosity results in a damping of the wave, we expect to have a frequency-dependent 
absorption of GWs in this frequency range. A spectral feature in the nHz range would 
be potentially interesting for PTAs. PTAs combine the fact that pulsars are very 
stable clocks with the fact that the time of arrival of the pulse also depends on the 
GW background between the pulsar and the Earth [TTJ. The use of a pulsar array 
allows to correlate the electromagnetic signals from different pulsars, thus eliminating 
variations in the arrival times that do not depend on the GW background. The 
dominant contribution to the GW background in the PTAs sensitivity range is that 
of supermassive black hole binaries following galaxy mergers [21] [22] [23] . This has 
an astrophyisical origin and would obviously not be affected by the decoupling of 
cosmological neutrinos. The dominant cosmological constribution to the GW signal 
in the nHz range is expected to be that of cosmic strings, that are topological 
defects left after symmetry-breaking phase transitions occurring in the early Universe. 

| In Ref. 1171 . the evolution of cosmological GWs in the presence of neutrino free-streaming was 
studied up to second order in perturbation theory. Moreover, GW propagation is also affected by 
entropy-generating events (like the e+e~ annihilation) occurring during the cosmological evolution; 
the effect on the GW spectrum has been studied in Ref. |18| . The authors of Ref. |18) also consider 
the effect of neutrino free streaming but assume an instantaneous neutrino decoupling. A similar 
analysis, including an accurate computation of the primordial spectrum as generated during inflation, 
and extended to higher frequencies, can be found in Ref. 19 . Finally, in Ref. 1201 . the role of a 
non-vanishing lepton asymmetry in determining the amount of absorption has been considered. (Note 
added in the arXiv version.) 
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Other, possibly important, contributions in this frequency range are those of cosmic 
superstrings (i.e., topological defects arising in string-theory inspired inflationary 
models) and of inflation-generated GWs (through the usual amplification of quantum 
zero-point fluctuations). For our purposes, one important difference between the 
cosmic (super)strings-generated GWs on one side, and the inflationary GWs on the 
other, is that the former are created by causal processes at horizon and subhorizon 
scales, while the latter are generated by the stretching of quantum fluctuations to 
super-horizon scales. In the following we will concentrate on GWs that were outside 
the horizon at some time before neutrino decoupling, so that our results rigorously 
applies to inflationary GWs but not necessarily to string-generated GWs. It should be 
noted that the inflationary signal is, in the nHz range, the weakest of the cosmological 
signals mentioned above, and is probably beyond the reach of the ongoing PTA 
projects. It could be however within the reach of future instruments [TTJ [53] like 
the Square Kilometer Array (SKA) |25j . 

The paper is organized as follows. In Sec. 2, we recall the basic equations 
describing the evolution of GWs in the presence of an imperfect fluid. In Sec. 3, 
we discuss the interaction between GWs and the neutrino fluid in the three regimes 
T > 1 MeV, T ~ 1 MeV and T < I MeV, corresponding to / > 1 nHz, / ~ 1 nHz and 
/ <C 1 nHz respectively. In Sec. 4, we present the results of the numerical integration 
of the Einstcin-Boltzmann system, show the modification in spectrum of cosmological 
GWs and briefly discuss the conditions that could lead to its detection. 



2. Basic equations 



We shall use, all throughout the paper, natural units in which c = h = ks = 1. 

Let us consider a gravitational wave, propagating on the background of a flat 
(k=0) Friedmann Universe. In synchronous gauge, the spatial components of the 
perturbed metric are written as [we use the ( — h H — (-) signature for the metric] 

9ij = a 2 (t)[8 ij + h ij ] (1) 

while the other components are left unperturbed: 500 = — 1 and goi = 0- We will 
consider only the transverse traceless part of With this restriction, only has 
two degrees of freedom, corresponding to the two polarizations of a GW. Here ait) 
is the cosmological scale factor, that evolves according to the background Friedmann 
equation: 



da\ 2 8ttG 



dt 



a 2 p (2) 



where p is the background density of the cosmological fluid (in general, we will use 
overbars to denote background quantities) . 

The Einstein equation for the evolution of a transverse traceless metric 
perturbation hij takes the form of a wave equation [21 [26l [27] : 

where is the anisotropic stress, i.e. the traceless part of the three dimensional 
energy-momentum tensor T' of the cosmological fluid, representing dissipative effects 
that are not present in a perfect fluid. It is then defined through the relation 
nj = Tj - <5)T*73 = T] - V8), where V is the total pressure of the fluid (including 



Signature of neutrino decoupling in the nHz region of the GW spectrum 



4 



possibly a small perturbation with respect to the background). We assume that the 
anisotropic stress is relevant only for one component of the cosmological fluid (the 
neutrinos), providing a density p v = f^p, while for the other components (e.g., the 
photons) it can be safely neglected. 

In the following, we will find more convenient to use the conformal time n, related 
to the synchronous time t by dt — a{rf)drj as our time variable. Moreover, since the 
perturbation equations are linear, it is very convenient to Fourier transform the spatial 
dependence of all the relevant quantities. Then the wave Eq. ((3]) rewrites as: 

hij + 2Uhij + k 2 h tj = WttGU^ , (4) 

where dots denote derivatives with respect to n, and H is the conformal Hubble 
constant % = a/ a. For the sake of simplicity, we use the same symbol for a given 
quantity and for its Fourier transform. 

In a perfect fluid, anisotropic stress is absent by definition (in fact, for a perfect 
fluid = V5j), so that Eq. (j3|) takes the homogeneous form: 

hij + 2Hhij + k 2 h zj = 0. (5) 

In the radiation-dominated era, a oc r\ and % = 1/r/. Then, considering the quantity 
Xij = kr/hij we find that 

Xij + k 2 Xij = 0. (6) 

This is a simple, undamped oscillator equation, that admits the solutions xij = sm kr] 

and Xif — cosfc?7, i.e.: 



hij = aij — 1- bij — — — . (7) 



sin(fc?7) cos(kn) 

In the limit of kn — > 0, i.e., for a wave that is far outside the horizon, we have that 

Then, the most general solution outside the horizon (during the radiation dominated 
era) is the superposition of a constant mode and of a mode that decays like I/77. This 
is actually true also in the presence of anisotropic stress, because dissipative effects 
cannot operate on scales larger than the horizon, so that II* has to vanish in the limit 
kn 0. 

The macroscopic properties of the cosmological fluid, including the anisotropic 
stress, can be derived by the phase space distribution of its particles. The phase space 
is described by three positions x % and by their three conjugate momenta Pi = mdxi/ds, 
and evolves according to the Boltzmann equation: 

where L = d/dn is called Liouville operator, and C is the collision operator accounting 
for changes in the distribution function due to collisions between particles. Using the 
geodesic equation, the total derivative of / can be rewritten as: 

df_ = dx^dj_ _ p df_ 

drj dr] dx» af} dP»' y 1 
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The link between the Einstein Eq. Q describing the evolution of the metric 
perturbation hy and the Boltzmann Eq. Q is given by the following expression 
for the energy-momentum tensor in terms of the distribution function: 

If P^P 
T£ = —j f(x\ Pj, t)-^ dP 1 dP 2 dP 3 , (11) 

where g is the determinant of the metric. Then, once the collision term is also specified, 
Eqs. (HJ), (HI) and (jlip are all that is needed, at least in principle, to follow the 
propagation of a GW. 

In general, the Einstein-Boltzmann system for the coupled evolution of hij and 
/ is an integro-differential system, because the source term in Eq. (J3J) is given by an 
integral over /. The system can be however reduced to an ordinary differential system 
in the case of massless particles. This is made by integrating out the dependence of / 
on energy and expanding its dependence on the direction of the particle momentum 
in Legendre polynomials. This method is widely applied in the numerical treatment 
of scalar perturbations 29 . The analogous equations for the tensor case have been 
derived in Ref. |28j . to which we refer the reader for further details; here we just quote 
the relevant equations. The coupled Einstein-Boltzmann system is equivalent to the 
following system of infinite differential equations in fc-space: 
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In the above equations, the .F^'s and the C^'s are related to the coefficients of 
the expansion in Legendre polynomials of / and C[f], respectively (see Ref. [28] for 
details). The anisotropic stress is given by: Ily = a 2 p l/ J r J ^/47r. 



3. Interaction of gravitational waves with neutrinos 

In this section we discuss how cosmological neutrinos affect the propagation of 
cosmological GWs. We restrict our attention to waves entering the horizon well 
before the time of matter-radiation equality, corresponding to a redshift z ~ 10 4 . 
This corresponds to waves with a present frequency v 3> 10~ 16 Hz. We recall that in 
the radiation-dominated era, a oc r\ and T~L = l/rj. 

The thermal evolution of neutrinos can be divided into two distinct regimes. 
In the early Universe, when the temperature was sufficiently high, neutrinos were 
kept in thermal equilibrium with the cosmological plasma by frequent reactions like 
vv <H> e + e~, ve O ve, etc. As the Universe expanded and cooled down, the interaction 
rate r„ for these reactions eventually became smaller than the expansion rate H and 
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the neutrinos decoupled from the cosmological plasma. From this point on, neutrinos 
can be considered as freely streaming through the Universe. The temperature T doc that 
separates the collisional from the non-collisional regime, defined as the temperature 
when the interaction rate is equal to the expansion rate, i.e. T„(Tdec) — H(T^ ec ), 
is found to be Td ec — 1 MeV, corresponding to z ~ 10 10 and t ~ a few seconds. A 
wave that enters the horizon at that time has a frequency / ~ 1 nHz. Higher frequency 
waves enter the horizon before neutrino decoupling, and viceversa. Then, the following 
three regimes for the evolution of a GW can be correspondingly identified: 

• When T>1 MeV, the neutrinos are tightly coupled to the cosmological plasma, 
that behaves like a single, perfect fluid. The anisotropic stress of neutrinos (as 
well as of the other components) is negligible and the GW evolves according to 
the homogeneous Eq. <j5j> . This regime is relevant for waves with / 3> 1 nHz. 

• When T ~ 1 MeV, the neutrinos are decoupling from the plasma. They cannot 
be considered as a perfect fluid because the mean free path of particles is getting 
large, so that the anisotropic stress cannot be neglected. However, the collisions 
have still to be taken into account. The GW evolution has to be calculated from 
the full system Eqs. (fT2| - (TT6|) with the C^'s modeling the collision processes. 
This regime is relevant for waves with / ~ 1 nHz. 

• When T -C 1 MeV the neutrinos are collisionless and behave as free particles. The 
anisotropic stress cannot be neglected: the GW evolution has to be calculated 
from the full system Eqs. (fl2 ]) -(fl6 ]l with C (i) = for all Ps. This regime is 
relevant for waves with /< 1 nHz. 

In the first, high-temperature regime, the metric perturbation simply evolves as 
hij — hij = Xij sm(kr))/kri (assuming the initial velocity vanishes), where A^- is the 
initial value of the perturbation. In other words, the GW oscillates with a time- 
dependent amplitude oc 1/krj. This can be seen as the energy loss of the wave due to 
the expansion of the Universe. 

The third, low-temperature regime was first studied in Ref. |14) . where it was 
realized that the anisotropic stress of neutrinos would lead to a partial absorption of 
the GW. In particular, for the standard value of the neutrino density, it was found 
that the amplitude of the wave when it is well inside the horizon is reduced to 80% of 
its value in the absence of stress. The intensity is reduced to (0.8) 2 = 64% of its value 
in the absence of stress. In Ref. [14] . the Boltzmann-Einstein system was treated in its 
integro-differential form. We re-analyzed the problem in Ref. [28] using the multipole 
formalism, and allowing for different values of the neutrino density. In general, one 
has that the behaviour of the wave well inside the horizon [kr] 3> 1) is of the kind: 

^=TO / in( y^ (fen » i), (17) 

where D is a damping factor, is a phase and /„ = p v jp is the neutrino fraction, 
i.e. the ratio of the neutrino density p v to the total density p. This is constant 
in the radiation-dominated era since both p„ and p scale like a -4 , and is equal to 
/„ = 0.4052 for the standard case of three neutrino families with a present temperature 
T® = (4/ll) 1 / 3 T° = 1.9 K. In the following we will neglect the phase shift <j>, since we 
are only interested in the behaviour of the wave when averaged over many periods. 
The damping factor D only depends on f v . The result of Ref. [14] can thus be 
written as D(f v = 0.4052) = 0.8. In Ref. [28] we have found the slightly larger value 
D{f v = 0.4052) = 0.88. We have also calculated the maximum amount of damping, 
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corresponding to -D(/„ = 1) = 0.75. It is worth stressing that the cosmological 
neutrino background has not been directly observed yet. Although strong deviations 
from the standard scenario with /„ = 0.4052 are unlikely, the possibility that /„ has a 
different value (due for example to the presence of additional free-streaming particles 
in the early Universe) should be taken into account. 

It is worth stressing, at this point, that the behaviour of the wave in either of the 
two extreme regimes / 3> 1 nHz and / < 1 nHz does not depend on the frequency of 
the wave nor on the values of the cosmological parameters (at least for waves entering 
the horizon before matter-radiation equality). This is due to the fact that, when 

= and a oc r\ the explicit k dependence in the Einstein-Bolztmann system can 
be eliminated by considering the time variable u = krj. Also, the convenience of using 
u is that u = krj ~ 1 corresponds to the time of horizon crossing. Moreover, by writing 
hu — \ijh(rj), it can be seen that h(ji) satisfies the same differential equations as h^ 
[upon a redefinition of the J*^'s and C^'s], so that we can always take without loss 
of generality \j = 1. Taking all this into account, we have that the GW evolution for 
kr\ 3> 1, in terms of Xij = krjhij = uhij is: 

( sin(u) (/» InHz) 

{ £>(/„) sm(u) (/ < 1 nHz) 

The fact that the (rescaled) amplitude of the wave is going from 1 for / > 1 nHz 
to D(f v ) < 1 for / < 1 nHz already points out to the fact that one should expect a 
frequency dependence in the region / ~ 1 nHz. In order to quantify this, one should 
solve the evolution equations for waves entering the horizon in the intermediate regime 
where neutrinos are neither a perfect fluid nor can be considered as freely streaming. 
The effect of collisions has to be taken into account. The exact computation of 
the collision term C[f] (and consequently of the C^'s) depends on the details of 
the interaction. However, a useful although rough approximation consists in writing 
C[f] = —Sf/r, where 5f = f — f^ is the deviation of the neutrino distribution 
function / from its thermal equilibrium value f^°\ and r is a characteristic time of the 
interactions that mantain the equilibrium. This form, albeit very simple, captures the 
main features that we expect from the actual collision term. When the characteristic 
time of the interactions is very large, t — > oo and C[f] —¥ 0, i.e. the fluid is collisionless. 
On the contrary, when the characteristic time is very small r — ¥ and C[f] is very 
large unless 5f — > too, i.e. the frequent collisions tend to mantain the equilibrium. 
This form for C[f] gives [3S]: 

Cj? = -3® It. (19) 
In our calculations, we take r to be the mean time between collisions, given 
by t = (^(crt;))^ 1 , where n is the number density of particles and (cry) is the 
thermally averaged cross section times velocity. The number density of particles in 
thermal equilibrium is n ~ T 3 , while an adequate value for the cross section for weak 
interaction processes is a ~ G^T 2 , where Gf is the Fermi constant. Then: 

t ^-gP^- (20) 

4. Results and discussion 



Equipped with the form (TT51) for the collision terms, with r given by Eq. (|2U|) . we 
have integrated the Boltzmann-Einstein system for 71 values of / = k/2n between 
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10~ 12 Hz and 10 _85 Hz ~ 3 nHz, with a constant logarithm spacing A log / = 0.05. 
We have assumed the standard value f v = 0.4052. For each value of /, we have defined 
a frequency-dependent damping factor Df as: 

D f = ^{{knKjY) (21) 

where the brackets (. . .) denote the average over many periods, taken when the wave 
is well inside the horizon (kr) 3> 1). It is clear that when h{j — sm(kr))/ki] we have 
Df = 1. The damping factor basically quantifies how much the amplitude of a GW 
produced in the early Universe would be reduced due to the neutrino anisotropic 
stress, with respect to propagation in a perfect fluid. We express our results in terms 
of the density parameter of GWs J7gw(/) = (dpow / dha /) / p c , where pgw is the 
energy density of GWs and p c is the critical density of the Universe. The quantity 
J7gw is used by theorists to quantify the intensity of the cosmological GW background 
predicted in a given scenario, and experimental sensitivities are also often quoted in 
terms of it. The present intensity is usually computed simply by rescaling the GW 
amplitude at the source by a factor 1/krj to take into account the redshift due to the 
expansion. This should be corrected to take into account also the anisotropic stress of 
neutrinos. Since pew °c h?j, the correction factor is simply D^. Thus, if we call f2cw 
the value in absence of stress, we have that Slew = D ?Ogw- I n Fig- fflwe show 
as a function of frequency. 



1.00 r-r 




0.70 u 1 . . — i — , , , 1 . . — i — ■ ■ ■ ■ i 1 J 

1 x 10" n 2x 10" 11 5x 10" n l x 10" 10 2x 10" 10 5 x 10" 10 1 x 10" 9 2 x 10" 9 

Frequency [Hertz] 

Figure 1. The correction factor D 2 j = QgwI&GW- The highest frequencies in 
the plot correspond to waves entering the horizon when the neutrinos are tightly 
coupled to the plasma, so that there is no dissipation. The lowest frequencies 
correspond to waves entering the horizon when the neutrinos are freely streaming, 
so that D 2 is constant and roughly equal to 0.75. 

The fact that the collisions introduce a frequency dependence on the amplitude 
of the GW background (in addition to any dependence that could be already present 
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at the source) is a more interesting case with respect to a frequency independent 
suppression, since the latter could not be disentagled from our ignorance on the 
amplitude of the original spectrum (that usually depends on the free parameters of the 
theory and thus cannot be determined a priori without any additional experimental 
input). On the other hand, a frequency-dependent effect can be disentangled by using 
this same dependence. To better quantify this, let us suppose that the GW spectrum 
at the source, has a featureless power law behaviour oc f a . After correcting for the 
expansion, the shape is still the same. Then the spectrum, once the effect of anisotropic 
stress has been taken into account, will be: 



n GW (f) = ADjr. 

The logarithmic slope of the spectrum in a given point is given 
d log VLqw I d log / , so that, defining the deviation Aa = a' — a, we have: 



(22) 



by a' 



Aa 



d log fi G w 



dlogD 



f 



dlogf dlogf 
In Fig. [5J we show Aa as a function of frequency, 
value of Aa = 0.15 at / ~ 0.1 nHz. 



(23) 



We find that it has a maximum 
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Figure 2. Modification Aa to the logarithmic slope of a featureless power-law 
spectrum f a at the source, introduced by the anisotropic stress. 



The nHz region of the spectrum of GWs can be probed by PTAs, combining 
the fact that pulsars are very stable clocks with the fact that the time of arrival of 
the pulse also depends on the GW background between the pulsar and the Earth 

Correlating the electromagnetic signals from different pulsars could, in principle, 
enable a positive detection of the GW background in the galaxy. PTAs are sensitive to 
GWs with frequency between roughly 1 and 100 nHz, with a maximum sensitivity in 
the region between 3 and 10 nHz. The lower bound comes from the fact that standard 



Signature of neutrino decoupling in the nHz region of the GW spectrum 



10 



pulsar timing techniques absorb any low-frequency signal, and so the time span of 
the data (currently ~ 30 years) gives a lower bound on the observable frequencies. 
The principal source of GWs in the region of maximum sensitivity are believed to be 
coalescing supermassive binary black- hole systems in the centre of merging galaxies. 
However, a signal of cosmological origin could also be present in that same frequency 
range, for example GWs produced during the inflationary era or from the oscillation 
of cosmic string loops. 

Can the effect of anisotropic stress leave an imprint that could be, in the future, 
detected using PTAs? This could be possible, at some conditions. First of all, 
obviously, a signal of cosmological origin has to be present in the nHz range and 
it has to be strong enough to be detectable by PTAs. As noted in the introduction, 
the oscillation of cosmic (super) strings loops could give a detectable signal in the 
nHz range. However, GWs generated in this way are produced at scales smaller or 
equal than the horizon, so that the analysis made here does not rigorously apply. 
Although it is sensible to expect that GWs with wavelengths comparable to the 
horizon would behave in a similar way to the super-horizon GWs considered here, 
we cannot assess the effect on the overall spectrum and we defer this generalization 
to a future work. This leaves the GWs produced during inflation by the amplification 
of quantum vacuum fluctuations, that arc too weak to be detected by ongoing PTA 
projects but could be within the reach of future experiments like SKA. In any case, 
an independent observation of the signal at larger frequencies, like those probed 
by interferometers, would also be useful to normalize the spectrum at the source. 
Secondly, the cosmological signal should be clearly separated from the astrophysical 
signal, like that produced by black hole binaries. The third point is the more 
problematic. We have shown that the effects of the damping are more evident at 
frequencies between 0.1 and 1 nHz. As explained above, the lower limit to the 
detectable frequency is given by the time span of the observations, currenly 30 years. 
In 70 years from now, with a total of 100 years of observations, the lower limit will be 
/ = l/(100years) ~ 0.3 nHz, where the intensity of the wave is reduced by only 5% 
and the deviation from a featureless power law is Act = 0.15. The maximum change of 
the slope is at / ~ 0.1 nHz, corresponding to 300 years of observations. It should also 
be noted that, even if we were willing to wait such a long time, it is not certain that it 
would actually increase the sensitivity, because of intrinsic low-frequency instabilities 
in the timing data. In conclusion, we think that the possibility to detect this effect 
are quite small, and that they will depend crucially on the existence of a sizeable 
inflationary-like GW background at the frequencies of interest and on the capability 
of estimating the logarithimic slope of the spectrum in the region just below 1 nHz. 
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